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ABSTRACT 


Two families of macroscopic adsorption models are formulated, based on fundamental 
principles of adsorption science and technology, that may be used for macroscopic (e.g., 
whole-building) contaminant dispersal analysis. The first family of adsorption models — the 
Equilibrium Adsorption (EA) Models — are based upon the simple requirement of equilibrium 
between adsorbent and room air. The second family — the Boundary Layer Diffusion 
Controlled Adsorption (BLDC) Models — add to the equilibrium requirement a boundary layer 
model for diffusion of the adsorbate from the room air to the adsorbent surface. Two members 
of each of these families are explicitly discussed, one based on the linear adsorption isotherm 
model and the other on the Langmuir model. The linear variants of each family are applied to 
model the adsorption dynamics of formaldehyde in gypsum wall board and compared to 
measured data. 
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1 
INTRODUCTION 


The dispersal of airborne contaminants and, thus, the quality of air in buildings is clearly 
dependent upon the nature of air movement into, out of, and within building systems and the 
character and location of contaminant sources. As a result, researchers in the field have focused 
on the complementary problems of modeling and measuring the airflow aspects of contaminant 
dispersal and the emission characteristics of contaminant sources. It is also well-recognized that 
some contaminants, such as Radon-222 gas and it's radioactive decay daughters and nitrogen 
dioxide, undergo chemical or radiochemical decomposition while being convectively transported 
by airflow and, therefore, modeling the dispersal of these contaminants demands the additional 
consideration of the chemical or radiochemical mass transformation processes involved. 


It is not well-recognized, but has certainly been suspected, that adsorption and desorption 
processes can also significantly affect the dispersal of contaminants in building airflow systems. 
The evidence to support this suspicion is, however, pervasive. The lingering odors of winter 
clothes stored during the summer in the attic with mothballs, the familiar smell of grandparents 
home that persists for month's, even years, after their home remains vacant, the perfume- 
soaked blouse, and the noxious odor of clothes that were exposed to tobacco smoke or gasoline 
fumes that persists for days until washed. Few would doubt that adsorption and desorption 
phenomena can and do play a role in the dispersal of airborne contaminants in buildings but, 
from the point of view of indoor air quality analysis, should we suspect these phenomena to be 
significant? Can adsorption and desorption mass transport processes significantly alter the 
dynamics of the dispersal of indoor air contaminants? Under what circumstances should one 
suspect the effect to be significant and how may one rationally model the dispersal of a 
contaminant when these circumstances occur? 


This report attempts to demonstrate that adsorption and desorption phenomena have not only 
the potential to significantly alter the dynamics of contaminant dispersal in buildings but may, in 
some cases, completely overwhelm the dynamics described by flow dispersal alone. Working 
from fundamental principles of the adsorption science literature, two families of 
adsorption/desorption models are formulated that provide a means to identify the circumstances 
under which sorption phenomena may come to be significant. It is shown that the simplest 
members of each of these families provide analytical models that appear to reliably predict 
indoor air quality dynamics when compared to measured data. 


Some fundamental principles of adsorption science are presented in the second section of 
this report. Two families of single zone adsorption dynamics models are formulated in the third 
section and from the simplest member of each family two classes of adsorption element 
equations are extracted that may be used for general multi-zone contaminant dispersal analysis. 
Strategies to implement these using existing multi-zone contaminant dispersal analysis programs 


are also presented. The adsorption models are then applied to two problems of modeling the 
adsorption and desorption dynamics of formaldehyde in gypsum wall board, in the fourth 
section, corresponding to experimental tests conducted by researchers at Oak Ridge National 
Labs and the National Institute of Standards and Technology. Finally, in the fifth section 
conclusions are draw and some directions for additional studies are suggested. 


2 
FUNDAMENTALS OF ADSORPTION SCIENCE & TECHNOLOGY 


Some fundamental ideas from the literature of adsorption science and technology will be 
reviewed here. Most of this literature is directed toward industrial applications of adsorption 
that involve the purposeful use of specially prepared adsorbents (e.g., activated carbon, 
zeolites, and synthetic resins) placed in adsorption devices (e.g., well-mixed batch reactors, 
well-mixed flow reactors, and fixed-bed reactors) that are carefully designed and controlled to 
achieve high rates of adsorption [11, 13, 14, 15, 16, 20]. The emphasis here, on the other 
hand, is directed toward modeling adsorption phenomena that is largely unintentional involving 
materials that tend to act as absorbents, but are not specially prepared to do so, that are placed in 
buildings for purposes other than that of adsorption separation. Given this emphasis we shall 
have to adapt ideas from the available literature and must limit our consideration to the simpler 
models of adsorption phenomena until additional study warrants the application (and provides 
the requisite data) of more complete adsorption models. An attempt will be made, however, to 
consider adsorption models that are sufficiently general to allow their application to modeling 
purposeful adsorption (e.g., indoor air purification) in building airflow systems. 


Adsorption involves the separation of a substance from one phase, indoor air in the present 
context, and the accumulation of that substance on the surface of another phase, building 
materials, here. The nature of adsorption will depend on the nature of the adsorbent, the nature 
of adsorbate, and the interaction that may occur between them. After Slejko [20], we may 
distinguish four generic types of adsorption — exchange, physical, chemical, and specific 
adsorption. Exchange adsorption involves the exchange of one ionic species — the ionic 
adsorbate — with other ionic species attached to the surface of adsorbent and, therefore, involves 
electrostatic attraction of the adsorbate to the adsorbent. In physical adsorption the adsorbate is 
bound to the adsorbent by relatively weak intermolecular Van Der Waals forces. Chemical 
adsorption or chemisorption results from chemical reactions between the adsorbate and the 
adsorbent and, therefore, involves usually very strong binding forces associated with the 
sharing of electrons between the adsorbate and adsorbent. Finally, the term specific adsorption 
is reserved for adsorption phenomena involving binding of molecular groups, without chemical 
transformation, to specific functional groups on the adsorbent (e.g., polar adsorbates binding to 
polarized sites on the surface of an adsorbent). Specific adsorption binding forces fall in the 
spectrum between those for physical adsorption and chemisorption. 


The nature of adsorption is directly linked to the strength of the forces that bind adsorbate to 
adsorbent with physical adsorption falling at the weak-force end of the spectrum and chemical 
adsorption at the strong-force end. With this in mind, it follows directly that: 

j 


a) adsorption is typically an exothermic process with the heat of adsorption associated with 
physical adsorption expected to be small relative to that associated with chemical adsorption, 


b) physical adsorption may be expected to be a relatively reversible process while chemical 
adsorption may not, and 


(c) chemical adsorption may be expected to be associated with the formation of a single layer 
of molecules of adsorbate on the surface of the adsorbent while in physical adsorption the 
possibility of multiple layers is more likely. 


Reversibility and the heat of adsorption of a given adsorbate/adsorbent pair thus provide a direct 
indication of the molecular mechanism responsible for the adsorption. 


Adsorption Equilibria 


Physical and chemical adsorption processes may be more or less reversible! and can be 
represented as: 


Adsorbate + Adsorbent | Adsorbate-Adsorbent Complex 


where the degree of reversibility depends, in part, on the thermodynamic state of the adsorbate- 
adsorbent system and the nature of bonding in the adsorbate-adsorbent complex. For physical 
adsorption the bonding is relatively weak and, thus, this reversible process may be thought, 
instead, to be a reversible exchange between two alternative phases of the adsorbate — the free or 
gas phase and the adsorbed phase — where the adsorbent may be thought to activate the phase 
transformation, not unlike a catalyst, rather than reacting to form a new adsorbate-adsorbent 


complex. Distinguishing the concentration of adsorbate near a surface, C , from that adsorbed 
on the surface, Cs, this reversible physical adsorption process may be represented as: 


adsorbent 
C Cs 


where, for our purposes, concentration will be expressed in terms of mass fraction (i.e., the 


concentration in the air phase, C , will be expressed in terms of mass adsorbate per mass of the 
gas phase (e.g., mass adsorbate/mass of air) while the concentration in the solid phase, Cs, 
will be expressed in terms of mass adsorbate per mass of solid adsorbent? ). 


1 For chemical adsorption reversibly may demand elevated temperatures. | 


2 The concentration of adsorbate adsorbed on the surface is sometimes expressed in terms of mass adsorbate per 
unit of specific surface area of the adsorbent or, in indoor air quality literature, in terms of mass adsorbate per unit 
of projected surface area of the adsorbent (e.g., wall surface area), Effective adsorbents typically have extremely 
large specific surface areas due to convoluted surfaces and porosity. The specific surface area of an adsorbent may 
be difficult to determine and may not be well correlated to projected surface area, hence, it is practically more 


Under isothermal conditions the rate at which adsorbate molecules (or ions) are adsorbed, in 
a Closed system, will eventually equal the rate at which they are desorbed, a state of equilibrium 
will be achieved, and the concentration in each phase will remain constant at their respective 


equilibrium values, Cg and Cse. The equilibrium surface concentration generally increases 
with the equilibrium gas phase concentration for a given isothermal condition. The variation of 


this concentration, Cse, with the gas equilibrium concentration, or for a given state of 
temperature, is commonly called an adsorption isotherm. Experimentally determined 
adsorption isotherms, typically represented graphically, may be approximated by one of several 
equilibrium models including the Linear Model, Langmuir Model, BET Model, and Freundlich 
Model; the first three having theoretical bases and the last being empirical. 


Linear Model: This adsorption isotherm is defined by a simple linear relationship: 


Cso = Kp Co (2.1) 


where Kp is a temperature-dependent coefficient known as the Henry adsorption 
constants or the partition coefficient [16, 20]. This single parameter (i.e., Kp) model 
may be expected to be appropriate for cases of physical adsorption of gases at very low 
concentrations (e.g., trace amounts) on homogeneous surfaces with ample sites available for 
adsorption. 


Langmuir Model: In the years 1916-1918 Langmuir developed an adsorption isotherm 
model that accounted for the fact that the sites available for adsorption are limited [16, 20]. 
The Langmuir model may be expressed, for our purposes, as: 


Cse = Eso Kt Ge (2.2) 
1 + Ki Ge 


where Cgpo is the surface concentration corresponding to complete coverage by a single 
monolayer (constant, for a given adsorbate-adsorbent system, and independent of 
temperature for the Langmuir assumption of a fixed number of sites available for adsorption) 
and K,_ is the temperature-dependent Langmuir adsorption coefficient. 


The Langmuir model approximates the linear model for very low concentrations of adsorbate 


in the gas (or liquid) phase (i.e., when K. Cr << 1) thus we should expect near equality 
between the partition coefficient and the product of the parameters of this model as: 


convenient to employ mass fraction. 

3 Henry's Law defines a linear, equilibrium relation between the vapor pressure of a slightly soluble gas and its 
concentration in liquid solution. By analogy to this gas solution process the linear model for adsorption is 
sometimes called Henry's Law. 


Kp = Cgo Ki (2.3) 


BET Model: Brunauer, Emmett, and Teller developed a model that accounted for the 
possibility of multilayer adsorption that may be expected to occur in physical adsorption 
[16]: 


Cso Keet (Ce/ ct. 


[1 -(Ce/o Ilr (Ceo, )+ Keerl*/o;. 


where KgeT is the BET constant related to the energy of adsorption and Cgo is the surface 
concentration corresponding to complete coverage by a single monolayer (constant for a 


Cse = (2.4) 


given adsorbate-adsorbent system) and one is the air phase concentration corresponding to 
saturated conditions. The BET model has been found to be reliable in the range of reduced 


concentrations of 0.05 < (oy Cla) < 0.35 [16] for many gas adsorbate-solid adsorbent 
systems. 


The BET model also approximates the linear model for very low concentrations of adsorbate 
in the gas (or liquid) phase: 


Cy Chal <<1 and Keer (Ce/ Aaa << 1 


thus we should expect near equality between the partition coefficient and the product of the 
parameters of this model as: 


Kp ~ (Cso Keet/o* (2.5) 


Freundlich Model: An alternative two parameter model that often fits experimental data well 
was promoted by Freundlich [20]: 


*(1/ 
Cse a Ke om m (2.6) 


where Kr and n are the constant parameters associated with the Freundlich model. 
From a somewhat simplified consideration of adsorption thermodynamics it may be shown 


that the temperature dependency of the adsorption coefficients used in the Linear, Langmuir, 
and BET models may be expected to be described by Arrhenius type relationships [16, 20]: 


Kp = Kp e-AH/RT (2.7a) 

KL = Ktpe-AH/RT C210) 

Kpet = Kpetoe(EteV/RT (2.7c) 
where, © 


a) Kpo, Kio, and Keto are constants, AH is the enthalpy change associated with adsorption 
(J/kmol), and E; is the average heat of adsorption of the first layer adsorbed which depend 
on the particular adsorbate-adsorbent system being considered, 


b) Eo is the heat of condensation of the adsorbate, and 


c) R is the universal gas constant (8314.4 J/kmol-°K), and T is absolute temperature (°K) at 
which a given adsorption isotherm is defined. 


A plot of the logarithm of experimental determined model constants log(Kp), log(KL), or 
log(Kpet) versus 1/T provides an indication of the correctness of the particular model (i.e., it 
should be a straight line) and provides a means to evaluate the Arrhenius constants. Borrazzo 
recently reported this to be so for some nine hydrocarbons commonly found in indoor air 
adsorbed on nylon, wool, and glass fibers [1]. 


3 
ZONE ADSORPTION DYNAMICS 


In this section we shall adapt some of the principles of adsorption science and technology to 
the problem of modeling adsorption dynamics in building zones. We begin by considering the 


adsorption dynamics of a single contaminant species a in a well-mixed building zone containing 
a volume of air with mass Mz and a quantity of adsorbent material with mass Ms. The 


(average) concentration of the contaminant species in the zone air will be identified as “C> and 
the average concentration in the adsorbent as “Cg , both expressed in terms of mass fraction 
(i.e., (mass &)/(mass air) and (mass «)/(mass adsorbent) respectively). 


Adsorbent Material 
“Cs, Ms 


a 
Building Zone = “Czw 
“Cz, mz 


Figure 3.1 Zone plus adsorbent control volume. 


If the mass flow rate of air into the zone is w and the species concentration flowing into the 
zone is “Cin, then the mass flow rates of species & into the zone and out of the zone, “Wj, and 


“Wout, are simply equal to the product of the respective concentrations and the air mass flow 
rate as: 


OWin = “Cin W (3.1) 
CWout = Cz Ww (3.2) 
In an analogous manner we may describe the species mass flow rate by adsorption and 
desorption transport between the room and the adsorbent material in terms of discrete or lumped 


species mass flow rates, “Wads and “Wes, that will, in general, depend on the zone and 


adsorbent concentrations as well as the physical characteristics of the adsorbate species Q, 
adsorbent material, and the detailed nature of the airflow in the room: 
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Wads feist Gz, "Cs, tee ) (523) 


OWdes = fides (“Cz, 2 Csitett ) (3.4) 


Isolating the control volume surrounding the adsorbent from the control volume surrounding the 
zone air, the adsorption and desorption mass flow rates may be diagrammatically represented as 
shown in Figure 3.2. 


Wout 


Fig. 3.2 Zone control volume and adsorbent control volume. 


Limiting consideration for the moment to completely reversible adsorption (i.e., allowing no 


irreversible destruction or transformation of the adsorbate species a in the adsorbent material or 
in the building zone) we may immediately write species mass balance relations for the adsorbent 
material, the building zone, and the systems as a whole as: 


Building Zone Mass Balance 


a - Ow. a ty d“Cz pa O. 
(“Wout - %Win) + ( Wads - “Wdes) hare ate Oy (3,5) 


Adsorbent Material Mass Balance 
d“Cs a 


(2Wdes - *Wads) + San Seis (3.6) 
Buildi Adsor. Material Mass B 

a a 
(2Wout - Win) + mz 2 47 mses = “G7 + “Gs (37) 


where we admit the possibility of generation of the contaminant species © in both the building 


zone, “Gz, and the adsorbent material, *Gs. 


4 


Irreversible destruction or transformation of the contaminant species in the adsorbent or in 
the building zone may be introduced by introducing appropriate transformation kinetics 
expressions for either or both of the generation terms above, “Gz and “Gs. Thus, for 
example, if the analyst wishes to model linear destruction kinetics in the adsorbent “Gg would 
be replaced by an expression of the form “Gg = (— mgks) “Cs, where ks is the first order 
rate constant associated with the destruction. To allow a focused consideration of (reversible) 
physical adsorption irreversible transformation will not be considered further. 


Equilibrium Adsorption 


To provide some insight into the dynamics of this coupled system it is useful to consider a 
limiting case where the adsorbed species concentration remains in equilibrium with the zone 
concentration (i.e., the species mass transport rate by adsorption or desorption is practically 
instantaneous). More specifically, if the equilibrium relationship is described by the Linear 
Model: 


“Cs =Kp“Cz Linear Model (3.8a) 
or: 
d“Cs d“Cz 
RaGs sei sdabz 
at at Ce 


the mass balance for the system as a whole, Equation 3.7 with Equations 3.1 and 3.2, 
simplifies to: 


C367 aie 
OFZ 206 
at 


w"Cz + (mz +Kpms) 


Linear Equilibrium Adsorption (3.9a) 


CE = (“Gz + “Gg + wCin) (3.9b) 


where all forcing terms have been collected on the right hand side to define the system 
excitation , “E. 


Alternatively, if the equilibrium relationship is described by the Lan gmuir Model: 


AC, = Csoki.%Cz 


1 + Ki°Cy Langmuir Model (3.10a) 


2 


Or: 


d"Cs _{_ %CsoK___| d*Cz (3.10b) 
dt (1 + K,°CzF dt 


the mass balance for the system as a whole simplifies to: 


a. a 
w%Crz + [me + eaese sae eee = ~E | Langmuir Equilib.Adsorption (3<41) 


(1 + Ki"CzF 


For both cases the adsorbent acts to increase the participating mass of the building zone from 
that of the mass of the air in the zone, Mz, to this quantity plus the mass of the adsorbent, ms, 
scaled by a factor. For Linear equilibrium, the effective participating mass is: 


Metfective = Mz + KpMs Linear Equilibrium Adsorption (312) 


For the Langmuir Model the scaling factor includes a nonlinear dependency on the zone 
concentration: 


Langmuir Equilibrium Adsorption (3.13) 


¥CaK 
Meffective = mz + (css |e 


(1 + K.%C2) 


but for situations where the zone concentration does not vary greatly, this factor will remain 
reasonably constant as well. Furthermore, for those situations where the zone concentration is 
small enough so that the KL“Cz term in the denominator may be neglected then this factor 
should be expected to equal the partition coefficient of the Linear Model (as discussed earlier). 


For a single well-mixed zone under conditions of constant flow the system time constant, T, 
(or its inverse, the nominal air change rate); 


t = (MZAy) ; W=constant Well-mixed Zone (3.14) 
provides a concise and convenient measure of the systems dispersal dynamics. For a well- 
mixed zone under conditions of constant flow with linear equilibrium adsorption dynamics 


added from Equation 3.9 it is seen that the system time constant, Tt, becomes: 


t = (MetfectiveAy) ; W=constant Well-mixed Zone with Equilib. Adsorption (3.15) 
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/// Example Application: Equilibrium Adsorption in a Room 

To fix these ideas and to evaluate the potential significance of the participating mass 
offered by real adsorbents in building construction it is useful, at this point, to consider a 
practical example. To this end, consider a room 5 x 8 m in plan with a 3 m floor to ceiling 
height, 25 mm thick nylon carpeting, and 13 mm thick gypsum board walls and ceiling, as 
sketched below. 


Gypsum Wallboard 
@ 10.0 kg/m? 


w =1.44x 10° g/hr 
(1.0 ACH) 


Nylon Fibers 
@ 2.5 kg/m? 


The approximate weight of the gypsum board may be estimated at 10 kg/m2 and that of 
the nylon fibers in the carpeting at 2.5 kg/m2 thus, given the room geometry, the total mass 
of gypsum board is approximately 7.8 x 10° g or: 


and for the nylon fibers in the carpet approximately 1.00 x 109 g or: 
Mnylon = 1.0 X 10° g 


At standard conditions of temperature and pressure the density of the air in the room is very 
nearly 1.2 kg/m? thus the mass of air in the room is some 1.44 x 105 g or: 


mz = 1.44x10°g 


Researchers at both the Oak Ridge National Laboratories and the National Institute of 
Standards and Technology report partition coefficients for the formaldehyde 
(HCHO)—gypsum wallboard system at 23°C and 50% RH to be approximately 5.5 g-air/g- 
gypboard4 [12, 18, 19]. If we assume equilibrium adsorption conditions prevail, then, 


4 ORNL researchers reported a "CH,O capacity" of 23 (m3-air/m2-gypboard) while Silberstein at NIST 


measured capacities of 46 and 47 (m3-air/m2-gypboard). Both measurements were based on mixed units of 
HCHO surface concentration expressed in units of (mg-HCHO/m2-gypboard) and air concentration expressed in 
units of (mg-HCHO/m3-air), but ORNL considered exposure of both sides of gypsum board specimens while 
Silberstein considered exposure on only a single side. When expressed in mass fraction units these effective 
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for this measured partition coefficient and the representative room considered here, the 
gypboard may be expected to contribute a participating mass of: 


MEMO-OWPPPAKE Mayy = 5.52 (g-ait/g-gyp) x 7.8 x 10° (g-gyp) 
= 43 x 10° (g-air) 


that effectively increases the room air mass, Mz, by a factor of 31. This effective increase in 
room air mass would result in an equal increase in the apparent time constant of the well- 
mixed zone, for conditions of constant airflow. Thus, for example, for a well-mixed zone 
without adsorption having a time constant of 1 hour (i.e., an air exchange rate of 1 air 
change per hour) the gypboard would effectively alter the time constant for HCHO dispersal 
to 31 hours. 


In another study [1] partition coefficients for nylon, wool, and glass fibers are reported 
for some nine different aromatic and chlorinated hydrocarbons, designated simply as H-C- 
Cl here. For nylon fibers these partition coefficients range from 2.4 to 16.8 x 10-4 g-air/g- 
nylon fibers>. If, again, we assume equilibrium adsorption conditions prevail, then, the 
nylon fibers in the carpet of the representative room being considered may be expected to 
contribute a participating mass of: 


eerVe"KE Maylon = (2.4 to 16.8) x 104(g-air/g-nylon) x 1.0 x 10° (g-nylon) 
= (2.4 to 16.8) x 10 (g-air) 


In this case the room air mass would be effectively increased only slightly by 0.01 to 
0.10%. /// 


These simple calculations based on the equilibrium adsorption assumption indicate that 
adsorption can dramatically alter the dynamics of contaminant dispersal in rooms when the 
product of the adsorbent mass and its partition coefficient are significant. Some first 
experimental studies of adsorption in rooms support this conclusion [12, 18, 19]. 


Boundary Layer Diffusion Controlled Adsorption 


In general, adsorption mass transport rates are determined by the individual transport rates 
of the diffusion processes and the surface adsorption kinetics that together account for the 


partition factors are, therefore, equivalent. Conversion from these units to the mass fraction units of (g-air/g- 
gypboard) was based on an air density of 1.2 kg/m, gypboard bulk density of 769 kg/m3, and a gypboard 
thickness of 13 mm. 

5 The reported partition coefficient of 0.2 to 1.4 (cm3-air/g-nylon) were converted to the mass fraction units of 


(g-air/g-gypboard) based on an air density of 1.2 kg/m. 
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multiple steps in the overall adsorption transport process. Throughout the adsorption literature, 
however, it is emphasized that adsorption dynamics, in most practical situations, is invariably 
controlled by diffusion transport - the adsorption kinetics step is practically instantaneous - 
although exceptions to this general rule exist© [8]. With this in mind, we will consider 
models for room adsorption dynamics that account for external film diffusion based on the 
assumptions that: 


a) equilibrium exists at the solid-gas interface at all times as defined by one of the 
equilibrium isotherms considered above (i.e., the kinetics of adsorption is practically 
instantaneous) and 


b) the system remains practically isothermal (i.e., the rate of heat generation due to 
adsorption is negligibly small and the heat transfer processes are sufficiently large to 
maintain isothermal conditions). 


Models for room adsorption dynamics that account for macropore diffusion and micropore 
diffusion will be considered in subsequent studies. 


Adsorption mass transport rates from the room air to solid (i.e., nonporous) adsorbents with 
reasonably smooth surfaces may be expected to be controlled by the rate at which a contaminant 
species diffuses from the bulk phase to the surface of the adsorbent. The diffusion mass 
transport rate will, in general, depend on the details of the contaminant concentration and 
velocity fields in the room. Details that are both difficult to measure and to determine 
analytically (e.g., by computational fluid dynamics) and, consequently, details that will not 
normally be known nor easily determined. For wall surfaces and other larger flat surfaces in 
rooms it is likely that any flow that exists will be directed more or less parallel to these surfaces 
forming a boundary layer over the surface of interest. Therefore, to obtain estimates of 
diffusion mass transport for these surfaces, we may reasonably apply some key results from 
boundary layer theory. 


The adsorbent surface may be thought to be separated from the room air by a film of air 
(i.e., the boundary layer) over which the contaminant concentration varies from a near-surface, 
air-phase concentration, aC, to the bulk air-phase concentration in the zone, “Cz, as illustrated 
in Figure 3.3. (Anticipating the formulation, we will assume that this near-surface, air phase 
concentration will remain in equillibrium with the adsorbed-phase concentration and, hence, use 
the "*" notation indicating an equilibrium concentration.) 


6 Ethier and Kamm consider cases (e.g. the adsorption of large molecules on surfaces of very small pores) 
where adsorption kinetics controls the overall adsorption dynamics. 
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Adsorbent Material 
Ms 


Surface 


Roughness 
OWin Wout 


Building Zone HODES 
ers Mz Thickness “Cz 


Fig 3.3 Boundary layer diffusion controlled adsorption. 


If mass transport to the adsorbent is due only to film diffusion (e.g., there is no bulk flow 
of contaminant through the adsorbent due to pressure gradients within or across the adsorbent) 
and the air-phase concentrations are not varying too rapidly, then the net mass transport rate 
from the bulk phase to the surface, “Ws, may be approximated using steady-state mass transfer 


relations from boundary layer theory 7 [21] of the general form: 


Cws = ("Wags - “Wdes) = Am p As(*Cz . C3) (3.16) 
where; 
hm isthe average film mass transfer coefficient acting over the adsorbent surface 
(length/(time) 
p is the film density of air, the average of the bulk and surface densities 


As __ is the projected surface area of the adsorbent 


A dimensionless form of the average film mass transport coefficient, the average Sherwood 
number, Sh, , may be correlated with the surface flow Reynolds number, Re, and the air 
phase Schmidt number, Sc, as: 


Sh, = Om’ { 0.664Re,"2Sc'% ; for Re, < 500,000 (3.17a) 
a-airp 
Shc = Om“ = 0.037 Re **Sc'* ; for Re, > 500,000 ; Sc~1.0 (3.17b) 
Q-airp 
where; 


7 The correlations used here are based on steady conditions of surface and bulk concentrations, stationary flow 
conditions parallel to the wall, and diffusion due to concentration gradients only, ignoring the so-called blowing 
effect (which must be expected to be small for mass fraction concentrations less than 0.2). 
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Rey =. UIs 
Vv 


U is the mean fluid velocity (parallel to the surface) outside of the boundary layer 
L is the length of the surface in the direction of the flow 


Sc = —Y 
a-air 
Vv is the kinematic viscosity of the air phase 


Gai) is the molecular diffusivity of the binary a-air system being considered 


For many gases in air the molecular diffusivity is practically constant with concentration 
(although it varies directly with temperature and inversely with pressure) with values, at room 
temperature and atmospheric pressure, ranging from 10-5 to 10-4 m2/s and the Schmidt number 
is typically close to 1.0 (e.g., 0.2 to 2.0). For a Schmidt number equal to 1.0 the average 
Sherwood number varies as plotted below. 


1.0E+5 

1.0E+4 
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Sherwood 
Sarnnee 1.0E+2 
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1.0E-1 

1.0E+0 1.0E+1 1.0&+2 1.0E+3 1.0E+4 1.0E+5 1.0E+6 1.0E+7 


Rec 


At relatively low flow conditions the Sherwood number approaches 1.0 as molecular diffusion, 
rather than convective transport, becomes dominant. At very low flow conditions, however, 
the average Sherwood number, from these correlations, goes to zero when, in fact, molecular 
diffusion should remain dominant, therefore, for our considerations here we will assume: 


Sh, = Smt =10 ; for Re. <5: Sc=1.0 (3.17c) 


Q.-alrp) 


With the expression for film diffusion transport, Equation 3.16, in hand we may formulate 


equations that describe the system adsorption dynamics by first demanding equilibrium at the 
adsorbent surface: 


“Cs =Kp“Cz Linear Model (3.18) 
or 


“Cs ms “CooK “Cz 


- Langmuir Model (gly 
1+K_%Cz 


and then substituting these equations and Equations 3.1, 3.2, and 3.16, into the basic mass 
balance relations for the system, Equations 3.5 and 3.6, to obtain, after simplification: 


a) for the Linear Model: 
ee oO oO 
(“Cow -“Ciaw) + hm p Ag*Cz ~©S) mt Oz _ ag, (3.20a) 
Kp at 
ae a oO 
rib A *Cz -“Cs| yy med LsieetaGs (3.20b) 
Kp at 


or, alternatively: 


Linear Boundary Layer Diffusion Controlled Adsorption 


— hmpA 
Sahn DAS a rem fc 


(3.20c) 
es 


— hm pA 
Simp Acer 


where, we have collected all forcing terms on the right hand side as the system excitation 
vector, {*E}: 


a an. 
ey | cai ag (3.204) 
AG. 
b) for the Langmuir Model: 
a a 
a ar. a Aan, _ Cs d’C7 ec 71 
( Czw Cinw) +m p Ad Cz K (Ceo -Cs) ~8G,) +Mz i Gz (32214) 
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— a 04 
Bhp Real Ge e Se ee Os) so (3.21b) 
ey { ~” Ki("Cgo —"Cs) : : 


or, alternatively: 


Langmuir Boundary Layer Diffusion Controlled Adsorption 


(The Langmuir Model Equations 3.21 are valid only up to the condition of complete coverage of 
the surface (i.e., (Cs9 —Cg) > 0) at which point Equations 3.21 become singular.) 


The Linear Model and Langmuir Model variants of the system adsorption dynamics 
equations considered so far may be represented in concise notation as: 


[“W]{"C} + jem} = {"E} (3.22a) 
where: 
a) for Equilibrium Adsorption: 
{FC} = °Cz (3.22b) 
[“W] = w (3.22c) 
[“M] = mz + Kpms Linear Model (3.22d) 
or 
[“M] = et eee als Langmuir Model (3.22e) 
(1 + KL*CzP 


b) for Boundary Layer Diffusion Controlled Adsorption 


aC, \ 


OC, | (S:225) 


co) | 
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hel aim , 0 
[“M] = Fae | (3.22g) 
[“W] = r Linear Model (3.22h) 

—hmp As en kes, AS 

Kp 

or 

w+hm p As — _ lim pAs _ 

KL(Cso -Cs) 
[“W] = on Langmuir Model (37221) 
Ts hmpAs 
Ae 2m PASI 
eee Ki(Cega Ca) 


Note that the Langmuir Model introduces nonlinearity in the mass matrix, [“M], of the 
Equilibrium Adsorption (EA) Model (i.e., with respect to zone concentration, Cz) and in the 
transport matrix, [“W], of the Boundary Layer Diffusion Controlled Adsorption (BLDC) 
Model (i.e., with respect to zone concentration, “Cs). Similar models can be directly 
developed for the BET and Freundlich equilibrium adsorption isotherms. If developed, these 
models will also reveal nonlinearities in the mass matrix for the EA Model and the transport 
matrix in the BLDC Model although the forms or nonlinearity will be different from those 
above. These alternate models may be considered in subsequent studies. 


The essential dynamic character of the (scalar) EA Model is contained in the system time 
constant, t, or equivalently the system eigenvalue, A = 1/t, which is equal to 
A = [*M]"[“W] = (WmMerfective). Analogously, the essential dynamic character of the 
coupled BLDC Model equations is contained in the coupled system eigenvalues that satisfy the 
characteristic equation of the system: 


det] (“MJ [°w] -A[q1]| = 0 (3.23) 


where, for the Langmuir Model this eigenvalue problem would be evaluated at each specific 
state of the adsorbent concentration, “Cs, given the nonlinearity of the mass transport matrix in 
this case. In general this eigenvalue problem has two solutions, {A1, A2}; the inverses of the 
system time constants, {71, T2}: 


[A] = 


0 r2 0 1/T2 


eA 


If we substitute the Linear Model, Equations 3.20, into 3.23 we obtain of the familiar form: 


an? +br+c¢=0 (3.25a) 
with: 
a= 1.0 (3.25b) 
pee (“hmpAs + w)Kpmg + “hmpAsmz (3.25c) 
Kpmsmz 
— 
_ _mpAsw (3.25d) 
Kpmsmz 
The two solutions of this quadratic relation, A = (-b + — 4ac Vb? = 4ac Von, for the variables 
here, yield the two eigenvalues of the Linear Model: 
, _C#hmpAs + w)Kpmg + MhmpAsmz 
2Kpmsm ; 
haste (3.26) 
4 i epees + w)Kpms + “hap ASA 4otmpAsw 
a2 Kpmsmz "eaten 


It is interesting to consider three limiting cases: 


a) No Diffusion: When diffusion from the building zone to the adsorbent surface drops to 


zero, “hmpAs = 0, the system eigenvalues correspond to the eigenvalue of the well-mixed 
zone without adsorption, by direct substitution in Equation 3.26, as expected: 


A1 0 aa 
i = a) 0 | ; when: “hmpAs = 0 (527) 
Ors 0 0 


b) No Infiltration: When the infiltration into the zone drops to zero, w = O, the system 
eigenvalues become, by direct substitution in Equation 3.26: 


pam Ms + we ; 
baie Kpmsmz 
= ; when: W = 0(3.28) 
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Note that as the first eigenvalue increases linearly with the diffusion mass transfer rate, 


“hmpAs, with the result that the system time constant, t; = 1/A1, approaches zero as the 
mass transfer rate increases. 


c) Equilibrium Adsorption: If the rate of diffusion, “hmpAs, is large relative to the 
infiltration mass flow rate, w, then the system eigenvalues approach those corresponding to 
the Linear Equilibrium Adsorption Model (i.e., from Equations 3.12 and 3.15): 


(el Ms + mz) 


0 
110 Kpmsgmz mn: 
un ;when:“hmpAs >>Ww = (3.29) 
0 Ao 0 ee ee Wis See 
Kpms + MZ 


Again the first eigenvalue increases linearly with diffusion mass transfer rate. 
The proof of Equation 3.29 is rather indirect. The quantities (“hp As “ w) in Equation 
3.26 may be replaced by “hmpAs for “hmpAs >> wand the resulting equation simplified 


to yield: 


- Oh pAs (Kp Ms + mz) 


rX 
‘ 2Kpmsm 
PMs Mz (3.30) 
,ahmpAs(Kpmg+mz) /, 4 wKkpmgmz 
2Kpmsmz AhnpAs (Kpms + mz 
By forming a Taylor's expansion of the square root term: 
wKpmsgmz 
1 - 4—— 
GhnapAs (Kpms + mzy 

og Su. " Ala wKpmgmz ae 
GhnpAs(Kpmg+ mz’ ‘2/\%hmpAg (Kpms + mzF 

=j1- pm VUSPING I Zee : when:"hmp As >> W (331) 


OhapAs (Kpms + mz 


substituting this result into Equation 3.30 and simplifying we obtain the results above. 
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/// Example Application: Boundary Layer Diffusion Controlled Adsorption in 

a Room 

To fix these ideas and to evaluate the significance of the diffusion transport reconsider 
the room studied earlier using the Adsorption Equilibrium Model now applying, instead, the 
BLDC Model. This time we will limit consideration to the adsorption of formaldehyde, 
HCHO, by gypsum board having the partition coefficient considered earlier, Kp = 5.5 g- 
air/g-nylon fibers. The area of the adsorbent in this case is Ag = 118 m2. Assuming air 
flows uniformly from the inlet to the outlet, a reasonable estimate of the effective length of 
the adsorbent surface would be L = 8 m and an estimate of the free stream velocity of the 
air flow would be U = 8 m/hr = 0.00222 m/s (i.e., given the flow rate of 1 air change per 
hour, the volumetric flow divided by the room cross sectional area). Thus the effective 


Reynolds number would be (with the kinematic viscosity of the airofv =~ 1.5 x 1 0° m2/s): 


RefsuWi =) (0.0022) (8) tims 
Vv 1.5x 10° 


For a Scmidt number of 1.0 we, thus, obtain the average Sherwood number for this 
problem, from Equation 3.17a, as: 


Sh = 0.664 Re,'7Sc'* = 0.664(1173)'7(1.0)'% = 22.7 


For a representative gas diffusivity of HCHO-ait—y = 1.4 x 105 m2/s we obtain an estimate 
of the film diffusion mass transfer coefficient from Equations 3.17 as: 


1.4x 10° 


= OOn7 = 3.97x 10° m/s 


he as shi (SD) 


With these values in hand, an infiltration rate of Ww = 40 g/s (i.e., based on 1.0 air change 
per hour and an air density of 1,200 g/m3), and the zone and adsorbent masses computed 
earlier: 
mz = 1.44x10°g ; ms = 7.8x10°¢ 
we may apply Equations 3.25 to compute the eigenvalues for the gypboard/room system: 
Ay = 3.17x10%s1 = 1.14hr! or t, = 0.876 hr 


Ap = 1.15x10%s1 = 4.14x10% hr! or % = 242 br 


revealing a bimodal response characterized by a very long time constant of 242 hr (i.e., 
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10.1 days) and a quick response somewhat smaller than the nominal system time constant 
without adsorption (i.e.,(™Zw) = 1.0 hr). This sort of response behavior has, in fact, 
been observed in experimental studies [12, 18, 19]; the results of these studies will be 
considered subsequently. 


The response of this system, and other systems governed by the BLDC Model, is 
particularly sensitive to boundary layer diffusion rate as characterized by the average 
Sherwood number. Plotted below, for this example case, is the variation of the system time 
constants for a variation of average Sherwood number: 
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It is seen that the first time constant, T,, diminishes from 1.0 hr to zero as the second 
time constant, T2, drops from very high values for an average Sherwood number near 
zero® to an asymtoptic value for large values. In the present case this asymtoptic value is 
30.8 hours and is equal to the time constant we obtained for this case using the Equilibrium 

8 Recall, by equation 2.26 the second eigenvalue approaches zero for a average mass transfer coefficient, or 


equivalently, an average Sherwood number of zero thus the time constant becomes unbounded for a zero average 
Sherwood number. 


Ply 


Adsorption Model as predicted by Equation 3.29. /I/ 


We may conclude, then, that, at one extreme, with a near-zero rate of diffusion mass 
transfer, the system will respond to changes of concentration in the zone as a simple well-mixed 
zone with a time constant nearly equal to(Mz/w). For nonzero rates of diffusion the system will 
rapidly respond to changes of concentration with a time constant somewhat less than the simple 
well-mixed zone time constant then, if the second mode is excited (e.g., if contaminant is 
generated or concentrated in the adsorbent) a very slowly changing response will be observed. 
Diffusion in smaller macropores and/or diffusion in micropores (e.g., having pore diameters 
approaching molecular sizes) may be expected to have, in essence, a similar effect to that 
produced by very low mass transfer rates in the BLDC Model and may be responsible for this 
form of observed behavior. At the other extreme, very high rates of film diffusion mass 
transfer lead to the simple behavior predicted by the EA Model as expected. 


Surface Smoothness 

The BLDC Model presumes diffusion from the zone air to a smooth surface. This begs the 
question: what is meant by smooth? In the present context, it would seem reasonable to define 
an adsorbent surface as being "smooth" (i.e., sufficiently smooth to apply the BLDC Model) if: 


a) the surface roughness is small relative to the thickness of the boundary layer and 


b) the structure of the surface roughness is coarser than that considered to be macroporous 
(i.e., the width of surface irregularities or indentations are much larger that than the mean 
free path of the diffusing contaminant species) so that the diffusion may be considered to be 
controlled by Fickian diffusion. 


The thickness of the boundary layer, 5, may be estimated as [4]: 


6 = (3.32) 


db 
Shi 


For the example considered above we obtain 6 = any = 0.29 m which is certainly large 
relative to the (likely) physical roughness of the gypboard 
of the BLDC Model in this case. 


surface thus justifying the application 


Adsorption Element Equations 
Contaminant dispersal analysis theory may be formulated using an element assembly 


approach wherein equations approximating the behavior of the building airflow system or 
subsystem, the system equations, are assembled from equations that describe contaminant 
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mass transport in discrete portions or elements of the system model. The details of this 
approach are presented elsewhere [3]. In this section we shall simply describe element 
equations that correspond to the adsorption models used above in the single zone case that may 
be used for multi-zone system models of arbitrary complexity. 


In general, for the dispersal of species «, contaminant dispersal element equations have the 
following form: 


are 
jewe}(2C} + [ame ET — fagey 3.33) 


where {“C°} is a vector of contaminant species a discrete concentration variables associated 
with the element "e", a subset of system discrete variables {*C}, and {“g®} is a vector of 
element derived contaminant species generation rates. The element matrices [“w®] and [%m®] 
are square transformation matrices known as the element species transport matrix and species 
mass matrices respectively. 


The element equations corresponding to the models presented above may be derived directly 
from Equations 3.22 by simply subtracting the single zone mass and airflow transport terms 
from the system equations: 


Equilibrium Adsorption Element: a single node element, say element "e", with an adsorbent 


mass of m8 associated with a well mixed zone "i" with: 


{*C°} = (cP! (3.34a) 

[“we] = [0] (3.34b) 

[“me] = Kpm8[1] Linear Model (3.34c) 

or 

[ome] = i oe Langmuir Model (3.34d) 
(eakpeoe) 

{ge} = \7G8) (3.34e) 


Boundary Layer Diffusion Controlled Adsorption: a two node element, say element "e", 


associated with a well mixed zone "i" and an adsorbent "j" with an adsorbent mass of m8: 


Gg 


“CP | 


; i = well-mixed zone node, j = adsorbent node (3.35a) 


ron = | 


rs | 


[ome] = a Pil (3.35b) 


[“we] = (he p Ag} Kp Linear Model (3.35c) 


[owe] = (he, p Ag Ki(“Cgo -*CF) Langmuir Model (3.35d) 


vd 
= 


age} = 0 (3.35e) 
val 


font 


where all variables are defined as earlier (those variables with a following superscript "e" are 
specific to the adsorption element "e"). Again, note that the Langmuir Models introduce 
nonlinearities, now at the element level, and, for the boundary layer diffusion controlled case, 
are defined for adsorbent concentrations less than “Co. 


Implementation Using Existing Multi-Zone Programs 

The Linear Models for both cdses above may be implemented using existing multizone 
contaminant dispersal analysis programs that provide only well-mixed zones and discrete flow 
paths. For equilibrium adsorption modeling the analyst need only increase the associated well- 


mixed zone mass by Kpmé to account for mass transport to the adsorbent. 


The boundary layer diffusion controlled case requires a slightly indirect tactic. By 
redefining the element concentration variables as: 


eae | eh | (3.36a) 


the element arrays may be rewritten as: 


Pie, 0 0 
[~“mé] = | oueaieme | (3.36b) 
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SS (Ri p A8)| A ; a | Linear Model (3.36c) 


These element equations have the form of a subassembly of a well-mixed zone linked to inflow 
and outflow flow elements, thus, these element equations may be added to existing multizone 
contaminant dispersal models by linking a pseudo-well-mixed zone to the actual well-mixed 
zone containing the adsorbent by two pseudo-flow paths (i.e., flow elements) having equal and 
opposite flows. By setting the pseudo zone mass to Kpmé and the pseudo-mass flows to 
(he, e) A8) the behavior of the adsorbent will be modeled although the computed concentration 
response in the pseudo-well-mixed zone must be scaled by Kp to obtain the species 
concentration in the adsorbent, “C/. 


Zo 


4 
APPLICATION & VALIDATION 


In this section we will consider the results of adsorption/desorption studies of formaldehyde 
(HCHO) on gypsum wall board reported by researchers at Oak Ridge National Labs and the 
National Institute of Standards and Technology. Before considering the details of these studies 
another study deserves mention. Cleary and Sonderagger, at Lawrence Berkeley Lab, 
developed a model to predict the adsorption and desorption of water in wood construction of 
residential attics [6, 7]. The approach taken was very similar to that developed above for the 
boundary layer controlled diffusion adsorption model. They combined boundary layer 
solutions (i.e., to model the diffusion from the bulk water concentrations to the surface of the 
wood construction) with an equilibrium relationship between surface and adsorbed water 
concentrations to form an adsorption dynamics model. The equilibrium relationship used was 
not one of the those presented above, instead the investigators used an empirical fit of a cubic 
polynomial to equilibrium data. In comparisons of modeled response to measured response of 
moisture levels in an attic this model proved extremely accurate. 


ORNL Study 

Researchers at the Oak Ridge National Labs investigated the adsorption and desorption 
characteristics of gypsum wall board specimens in a small test chamber and demonstrated that 
an empirical three-parameter exponential response model could be fitted to measured response 
data with reasonable accuracy. A diagram of the test set up is shown below: 


Reference Gypboard 
hamber hamber 
Vol = 0.20 m3 Vol = 0.20 m3 


23°C, 50% RH 23°C, 50% RH 


Airflow: 
5.15 L/min 


Airflow: 
5.15 L/min 


Cin 
Fig. 4.1 ORNL Dual chamber test facility. 


In this investigation six specimens of gypsum board (0.25m x 0.24m) were placed in one of 
a pair of dual test chambers (volumes of 0.20 m3) and air, at 23° C and 50% RH, containing a 


controlled concentration of formaldehyde, EG was passed through this chamber and an 
identical reference chamber at constant and identical rates. The exit formaldehyde 
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concentrations from each of the chambers, “Cgyp and “Cref, were measured at periodical 
intervals and the ratio, R, of these concentrations, adjusted for initial conditions and normalized, 


were reported. Finally, three-parameter (i.e., A, Z, and t) exponential relationships of the 
following forms: 


R 


Z + (A-Z)(1 -e**) for adsorption (4.1) 


R = (A-Z)et* for desorption (4.2) 
were fitted to measured response data. It must be emphasized, to use the words of these 
researchers, that "Although exponential sorption and desorption processes may represent the 
physical behavior of the experimental system, the precise form of the models is empirically 
derived." 


The gypsum wall board chamber was modeled using the Linear variants of the Equilibrium 
Adsorption (EA) and Boundary Layer Diffusion Controlled (BLDC) models, Equations 3.9 and 
3.20 proposed above, and the reference chamber was modeled as a well-mixed zone, as 
diagrammed below: 


C ref 


Coyp-EZ 


Bey SS, 
bier SS 
Bory Ss 
os. ~ 

es 


Met = 20,1 06 


oo 
= . * 
S 6S a 
ea be 
hes ‘ 
3 
= . 
~ >, 5 


Ww = 371 g/h 
Reference Chamber Gypboard Chamber Gypboard Chamber 
Well-Mixed Model EA Model BLDC Model 


Fig. 4.2 Idealizations of ORNL test chambers. 


(For these and the following computational studies the density of the gypsum wall board was 
estimated to be 769,000 g/m3, the molecular diffusivity of formaldehyde in air at 23°C was 
estimated to be 1.39 x 10-5 m2/s using the Chapman-Enskog formula [5, 17], and the kinematic 
viscosity of air was estimated to be 1.53 x 10-5 m2/s.) 


The response of these modeled chambers was computed, using CONTAM87 [2], for a 
twenty eight day adsorption test followed by a four day desorption test. The results of these 
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analyses are compared below to the ORNL empirical model along with an indication of its 
uncertainty (i.e., the + err and - err curves): 


=> EA Model 
sx DC Model 


== ORNL 


— ORNL +err 
— ORNL -err 


0 2 4 6 8 10-7 12.144 16.4 18 |) 204 eta o meee 
Days 


Fig. 4.3 Comparison of modeled adsorption response and ORNL empirical model. 


«« DC Model 


= ORNL +err 
— ORNL-err 


Fig. 4.4 Comparison of modeled desorption response and ORNL empirical model. 
It is seen that the initial response appears to be well-captured by the EA model and the longer 


term response is better captured by the BLDC Model. An examination of the ORNL raw data 
(see Figure 3 of [12]) reveals greater scatter, especially in the bend of these curves, than the 
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error bounds shown above indicate, thus a response falling between the EA and BLDC Models 
would be considered practically accurate. 


The BLDC model is sensitive to the rate of diffusion mass transport that is correlated to an 
estimate of the flow Reynolds number. In the present case, the flow Reynolds number was 
based on an estimate of the mean flow velocity in the chamber equal to the volumetric flow rate 
divided by the square of the cube root of the chamber volume (i.e., an estimate of the chamber 
cross-sectional area) — more detailed information of the chamber geometry could not be 
obtained. This resulted in an estimate of Re = 4.25, an extremely low value; it is certain that the 
obstruction of the test specimens would have resulted in a larger estimate. The effect of higher 
Reynolds numbers (1.e., greater diffusion rates) would bring the BLDC Model curve closer to 
the EA Model curve in the early response and, thus improve the accuracy of the modeling. 


The desorption response falls between the EA and BLDC model predictions. Again, the 
accuracy of the estimate of the flow Reynolds number is critical; a larger and, presumably more 
reasonable, Reynolds number would have improved the BLDC modeled results by "pulling" the 
BLDC curve closer to the EA curve. 


The adsorption dynamics models presented in this report depend critically on the equilibrium 
adsorption isotherm model employed, in the present case, on the partition coefficient Kp. A 
partition coefficient of 5.5 (g-air/g-gypboard) was used for these studies and the NIST studies 
below. This is a value that was determined by both ORNL and NIST investigators using a 
rather intuitive method based on integrating the desorption response data that corresponds, in 
essence, to an integral form of the EA Model equations, Equations 3.9, or: 


to to 
| “C7 dt + (mz +Kpms)A%Cz -| OF dt (4.3a) 
t 


1 ty 


where one integrates response and excitation values over an arbitrary time interval, ty ... ta, and 
evaluates the change of concentration over the same interval; 


A"Cz = °Cz(te)— *Cz(ts) (4.3b) 
Solving this equation for the partition coefficient: 


to to 
| OF dt— “| OMe 
t 


ti 1 MZ 
Kowa he ee (4.4) 
ms A°Cz Ms 


defines the method used by these investigators to determine Kp (although both groups neglected 


of) 


to consider the second term on the right hand side, (mz/Ms), which was practically insignificant 
in both studies). 


This determination of the partition coefficient implicitly assumes linear equilibrium 
adsorption behavior and, therefore, stacks the deck in favor of the linear equilibrium adsorption 
model. Silberstein, in the NIST study, also, however, used an alternate strategy to determine 
the partition coefficient and obtained practically identical results thereby providing some 
validation of the value used here. It would be best to determine this equilibrium adsorption 
parameter or, better, the equilibrium adsorption characteristics using one of several available 
independent means to do so [16]. 


NIST Study 

In tests to study formaldehyde levels in a test house due to emissions from a variety of 
pressed-wood products, researchers at the National Institute of Standards and Technology 
found that it took days to reach equilibrium concentrations levels and suspected these long 
delays were due to adsorption and desorption of HCHO in gypsum wall board present in the 
test house [9]. To investigate this phenomena more directly these researchers conducted a 
separate investigation to study the adsorption and subsequent desorption of HCHO in a 
"medium-sized" dynamic measuring chamber, formulated a semi-empirical model of the 
dynamics of the sorption processes, and numerically extracted values for the model parameters 
from measured data [18, 19]. In distinction to the ORNL study, the environmental chamber 
employed was nearly room size, the gypsum board was attached to the walls of the chamber, as 
it would be in building construction, exposing: a single side, and the proposed model was 
based, qualitatively, on the likely physical mechanisms responsible for adsorption. This model 
was not, however, directly formulated from fundamental principles of adsorption science and, 
therefore, must be considered a semi-empirical model. 


In the NIST study a 1.2 m by 2.4 m specimen of 0.013 m thick gypsum wall board was 
secured to a wall of a "medium size" test chamber with inside dimensions of 1.22 m x 2.44 m x 
0.61 m. This test chamber was placed in a closed environment maintained at 23° C and 50% 
RH, a constant fresh air flow rate through the chamber of 1.01 air changes per hour was 
maintained, and, for adsorption studies, formaldehyde gas was generated within the chamber at 
a constant rate. 


The NIST chamber was modeled using the Linear variants of the Equilibrium Adsorption 
EA and BLDC models estimating the density of the gypsum wall board, the molecular 
diffusivity of the formaldehyde, and the kinematic viscosity of the air, as above, and employing 
the partition coefficient of 5.5 (g-air/g-gypboard). Diagrams of these models are shown below 
in Figure 4.5. The response of these models to an eleven day adsorption test, with 
formaldehyde generated within the chamber at a constant rate of 0.24 mg/h, followed by a 10 
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day desorption test, without the generation of formaldehyde, were computed using 
CONTAMS87. The results of these computational analyses are compared below to the measured 
data, Figure 4.6. 


Regrettably, there was some confusion regarding the actual formaldehyde generation rate 
and infiltration rate used during the test. Test data was obtained from two draft reports of the 
tests. In one of these reports an infiltration rate of 1.01 ACH was reported but the uncertainty 
associated with this value was great. The formaldehyde generation rate was reported to be. 
0.313 mg/hr for the adsorption tests yet hypothetical steady state concentration levels for an 
empty chamber were reported that indicated the actual generation rate was 0.24 mg/h; this latter 
value was used in the computations. Computations were based on an air exchange rate of 1.01 
ACH and a generation rate of 0.24 mg/h. (Silberstein suggested, in a personal communication, 
that an air exchange rate closer to 1.3 ACH and the generation rate of 0.313 mg/h might be more 
appropriate. Since the ratios of these values are nearly identical to the values used the computed 
results will be very similar.) 


As in the ORNL case the measured response very nearly falls between the EA Model and the 
BLDC Model responses. Again the BLDC response was especially sensitive to the estimation 
of the bulk mass transfer coefficient that was correlated with the flow Reynolds number. In this 
case the Reynolds number was estimated to be 110, a relatively low value. On Figure 4.6 is 
plotted the BLDC Model response for a ten fold increase in Reynolds number (i.e., a Y10 
increase in diffusion mass transfer rate) and it may be seen that the modeled response in this 
case more closely follows the measured response. It appears that a somewhat higher Reynolds 
number would have led to excellent agreement, but, at this point in time we can not provide a 
reliable means to estimate the flow Reynolds number accurately a priori and, therefore, the 
difficulty of estimating a representative Reynolds number must be considered a primary 
limitation of the BLDC Model approach. 
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Fig. 4.5 Idealization of the NIST test chamber. 
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Fig. 4.6 Comparison of modeled response and NIST measured data 


36 


s] 
CONCLUSION 


Two families of macroscopic adsorption models have been formulated, based on 
fundamental principles of adsorption science and technology, that may be used for the purposes 
of indoor air quality analysis. The first family of adsorption models — the Equilibrium 
Adsorption (EA) Models — are based upon assuming zone air concentrations and adsorbent 
concentrations remain in equilibrium at all times where the equilibrium condition may be defined 
by any number of equilibrium isotherms. Two members of this family are explicitly discussed, 
the first based on the linear isotherm model and the second on the Langmuir isotherm model. 
The second family of adsorption models — the Boundary Layer Diffusion Controlled 
Adsorption (BLDC) Models — are based on the use of boundary layer theory to model the rate 
of mass transfer from the bulk phase through the boundary layer to the solid adsorbent phase 
and the assumption that equilibrium exists at the solid-gas interface. Again this equilibrium may 
be modeled using any one of several adsorption isotherm models. Two members of this model 
are explicitly discussed, one based on the linear adsorption isotherm model and the other on the 
Langmuir model, as before. 


An analysis of the linear variants of the EA and BLDC models, presented in this report, 
demonstrated that the EA model and the simple well-mixed zone model (1.e., without adsorbent) 
are simply limiting cases of the BLDC model. This analysis provides formal proof of the 
physical consistency of the BLDC approach. 


It is shown that the linear variants of both the EA and BLDC models lead directly to 
contaminant dispersal sorption elements that may be used directly with existing contaminant 
dispersal analysis programs such as the CONTAM family of programs. These models have the 
attractive feature that they depend on a single parameter, the partition coefficient, that may be 
measured using one of several well-documented methods [16]. They must be expected, 
however, to be limited to the modeling of physical adsorption of trace contaminants because the 
linear adsorption isotherm upon which they are based is generally only valid for very low 
adsorbate concentrations. The BLDC captures the biphasic [18, 19] or bimodal nature that has 
been observed in adsorption and desorption tests but it relies on an estimation of the nature of 
airflow in the zone, as characterized by the flow Reynolds number, which, for room 
geometries, is a very uncertain problem. Diffusion through a boundary layer must be expected, 
however, to be an important step in practically all cases of adsorption and, short of complete 
solution of the microscopic equations governing the flow in a zone, the boundary layer 
solutions to this problem are the best available. 


Applications of the linear variants of the EA and BLDC models to formaldehyde adsorption 


and desorption tests were presented to provide some first validations of these simple models. 
The agreement between measured and modeled response was excellent. These practical 


3) 


applications and the example applications considered to clarify the development of the models 
clearly indicate the significant, in fact, overwhelming affect that simple physical adsorption and 
desorption transport processes may have on the dispersal of contaminants in buildings. 
Furthermore, physical adsorption must be expected to be pervasive depending not so much on 
special characteristics of adsorbate-adsorbent systems but simply on the availability of sufficient 
adsorbent surface. To date, sorption processes have been largely ignored in indoor air quality 
analysis and, importantly, in the closely related field of tracer gas analysis of building airflows. 
Both areas of analysis should be reevaluated in light of the implications of the results reported 
here. 


The success of the work reported by Cleary and his colleagues [6, 7] in modeling moisture 
variation in wood construction supports the hope that the nonlinear variants (e.g., Langmuir and 
BET adsorption isotherms) of the EA and BLDC models will prove effective in modeling 
physical sorption mass transport when adsorbate concentrations vary over wide ranges and may 
not be at trace levels (e.g., moisture variations in buildings in general), but this will have to be 
investigated. 


Not presented in this report are physical adsorption and desorption models based on existing 
knowledge of macro and microporous diffusion processes. It is clear, however, that 
macroscopic models could be readily developed based on finite element approximations to 
available governing differential equations for these classes of diffusion (e.g., the so-called 
Glueckauf Model for the diffusive transport of an adsorbate down a cylindrical pore in a 
macroporous adsorbent [10]). The numerical studies reported clearly indicate the potential 
importance of diffusion rate limitations — diffusion resistance can lead to extremely long time 
constants in adsorption dynamics. It.is entirely possible that the exponentially-diminishing 
"source" models reported by some investigators and the apparent irreversible, but small, 
conversion of formaldehyde in gypsum wall board are both due to more deeply adsorbed (e.g., 
within macro or microporous structures) constituents. It seems likely that a consideration of 
these more rate limiting diffusion steps in the adsorption process may lead to some rationally 
based physical "source" models. 
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